h = clf;
hold on
bcd = importdata('d2_bcd.0.1.log');
ms =  importdata('d2_ms.0.1.10r.log');
ms100 =  importdata('d2_ms.0.1.100r.log');

h_bcd = decimark(bcd(:,3), bcd(:,2), 5, 's', {'color','r','LineWidth', ...
                    2, 'markersize', 10});
h_ms = decimark(ms(:,4), ms(:,2), 5, 'd', {'color','b','LineWidth', ...
                    2, 'markersize', 10});
h_ms100 = decimark(ms100(:,4), ms100(:,2), 5, 'o', {'color','k','LineWidth', ...
                    2, 'markersize', 10});
%h_ms_z = decimark(ms(:,4), ms(:,3), 5, 'd', {'color','b','LineWidth', ...
%                    2});

legend([h_bcd h_ms h_ms100], 'BCD \alpha=0.1', ['Master/Slave \alpha=0.1 ' ...
                    '\rho=10'], 'Master/Slave \alpha=0.1 \rho=100');
xlabel('Time (seconds)');
ylabel('||A-WH||_F');
title('Regularized NMF: BCD vs M/S size(A) = 10000 \times 2000 k = 50');

print(h, '-depsc', '-r300', 'bcd_vs_ms_10000_50_2000.eps');